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Abstract. We give an elementary and rigorous proof of the Thomae type formula for the 
2 i singular curves fi N = IljLi(-^ — ^'2j) N _1 IIJLoO^ — -^2j+i)- To derive the Thomae formula we 

use the traditional variational method which goes back to Riemann, Thomae and Fuchs. An 
important step of the proof is the use of the Szego kernel computed explicitly in algebraic form 
, for non-singular l/./V-periods. The proof inherits principal points of Nakayashiki's proof 1311 . 

• obtained for non-singular Zjv curves. 

o . 

r ~" 1 ' 1. Introduction 

C | The original Thomae formula 136j expresses the zero values of the Riemann f?-functions with half 

integer characteristics as functions of the branch points of the hyperelliptic curve /i 2 = Yii=i 
that is 



X 

5-H 



o 

S : (1-1) %](0;n) 8 = (^A T ) 4 [](A 4fc - A 4i ) 2 (A Jfc - A,,) 2 , 

■ k<l 

1 c*j where e is a non-singular even half-period corresponding to the partition of the branch points 
Oh' {1, • • • , 2m} = {ii < • • • < i m } U {ji < ■ ■ ■ < j m }- The Riemann period matrix II and the matrix 
^ (— i of a-periods A = (§ a X 3 / n)i<i,j<m—i are computed in a canonical homology basis (a,f3). 

The modern interest in the Thomae formulae was initially stimulated by the finite-gap inte- 

Sgration of KdV and KP type equations where, according to the Dubrovin-Novikov program of 
" effectivization of the finite-gap solutions" , it was necessary to express winding vectors in the Its- 
. Matveev formula [2D] in terms of f?-constants (see e.g. ^01 EDI)- In this context Thomae formulae 

were recently considered in |12| . 

Thomae formulae were used to describe the moduli space of hyperelliptic curves with level 
two structure in terms of theta constants |30j . They give a generalization of the A function of an 
elliptic curve |30j (Umemura's appendix) and |13j . In |28j they were combined with the arithmetic- 
geometric mean to obtain criteria for the reducibility of ultraelliptic Jacobians. Thomae formulae 
also appear in a wider context of modern research on integrable systems, in particular in the 
Riemann-Hilbert problem and the associated Schlesinger equations (see, e. g-,[2I]). In this set 
of problems the Thomae formulae were used to give explicit expression for the r-function of the 
Schlesinger equation associated with hyperelliptic curve |21| . general curve with simple branch 
points |27| and non simple branch points Thomae formulas are also used to compute action 
variables in conditions of complete integrability (e.g. 9 ). In the conformal field theory Thomae 
formulae became relevant after the work of Knizhnik |25j who expressed correlation functions of 
the multi-loop string amplitude in terms of ^-constants of a Riemann surface with Zn symmetry. 
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A solution of Knizhnik-Zamolodchikov equations in terms of ^-constants was also obtained in |35| 
on the basis of Thomae formulae. 

After the classical Thomae paper the derivation of Thomae formulae in hyperelliptic case 
was given in many places: Fuchs [TJI, Bolza 0], Fay ^3], Mumford (3U]. But only recently 
Bershadsky and Radul 5, BR88 discovered a generalization of Thomae formula for Zn curve 
[i N = Yli—ii^ ~~ A») an d gave a heuristic proof of them on the basis of path integral formulation of 
the conformal field theory. Afterwards Nakayashiki developed a rigourous derivation of these 
formulae in the frame of the classical methods. 

In this paper we derive the Thomae formula for the singular curve 

m m 

(1.2) V N = Y[(^-hj) N - 1 l[^~X 2l+1 ), N>1 NeN, 

which has Zn symmetry and it will be called singular Zn curves. 

The modular properties of certain family of curves (|1.2|l were investigated by Burhardt [Jj], 
Hutchinson J^j, more recently by Shiga |35 and Koike Diez [7] and others. This curve (|1.2[1 
appeared in Zverovich |40j as the curve related to a solvable Riemann-Hilbert problem with quasi- 
permutation monodromy matrices. Explicit solution of this problem was derived in our recent 
paper (IT] which stimulated the present investigation. 

Our proof of Thomae formulae goes up to the original Thomae paper |36| and inherits principal 
steps of the Nakayashiki's proof 31) for non singular Zn curves. The non singular Zn curves are 
invariant under the group generated by the permutations of branch points, while the singular Zn 
curves are not. This is the main difference with the case treated by Nakayashiki. 

The main steps of our proof are the explicit algebraic formula for the Szego kernel associated 
with non-singular 1 /TV-period and the explicit formula for the canonical bimeromorphic differential 
kernel. The above kernel functions can be realized in ^-functional form, and in algebraic forms. In 
the former case the theory was developed by Fay 14- , in the latter case the theory has originated 
in the work of Klein \2'6i I24| and later it was developed by Baker and also ^BLEU and for 

nonsingular Zn curves by Nakayashiki |31| . The comparison of these representations leads to a 
certain relation which is the important point of the proof. We also use the notion of variation of 
branch points by means of Rauch formulae |33 | . |27 | . 

The paper is organized as follows. In Section 2 we are fixing the notations and remind basic 
notions from the theory of algebraic curves which we need for the foregoing development. In the 
Section 3 we describe characteristics of Abelian images of branch points. We introduce in this 
section families of non-special divisors supported on branch points and describe corresponding 
1/iV-periods. We describe canonical bimeromorphic differential and Szego kernels in the Section 
4 and develop their expansions. In the Section 5 we briefly discuss Kleinian formulation of the 
kernel forms. All these results are summarized in the Section 6 where the formulation and proof 
of the Thomae formula is given. We briefly discuss further perspectives in the last Section. 

2. The curve, its homologies, differentials, ^-functions and variations 
Consider the Riemann surface C of the curve 

(2.1) ^ N =p(X) q (X) N -\ 

m m 

(2.2) p(A) = n(A-A afc+ i), q(X) = l[(X-X 2k ). 

k=0 k=l 

The curve (|2.1|1 has branch points Pk = (A&, 0), k = 1, . . . , 2m + 1 and P2 m +2 — Poo = (oo, oo) 
and singularities at the points P 2 k , k — 1, . . . , m + 1. These singularities can be easily resolved [2H| 
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to give rise to a compact Riemann surface which we still denote by Cn,™,- The genus g of Cn,™ is 
equal to (N — l)m. 

A, (which we still denote by A) defines Cn.™ as a N— sheeted covering 



The projection (A, y) 
of the complex plane C 
non-branch point A G C 
by the action J : (A, y) - 



branched over the points P k , k — 1, . . . 2m + 2. The pre- image of a 



consists of N points. The iV-cyclic automorphism J of CM,m is given 
(A, py), where p is the A^-primitive root of unity, namely p = e 2 ^ . For 
P in a neighbourhood Ur of the point R £ Cjv,mj a local coordinate x(P), with a;(i2) = 0, is the 
function defined by 



(2.3) 



A(P) — A(.R), if i? is an ordinary point, 
X (P)= \ ty\(P)-X(R), if R = P k , fc = l,...,2m+l, 

if /\- 



The canonical homology basis, (oti, . . . , a(j\r-i) m ; /3i, ■ ■ ■ 
shown in the Figure Q Namely the cycles a.k+sm, k = 1, . 



P(N~i)m) € H(C,Z) of the curve C is 
. . , m lie on the (s + 1)— th sheet, s = 




first sheet 
s-th sheet 
N-th sheet 



Figure 1. The homology basis. 



0, . . . ,N — 2 and encircle anti-clockwise the cuts U^ =1 (A2j-i, ^2j)- The cycles fik+sm, k= 1, . . . , m, 
s = 0, . . . , N — 2, emerges on the (s + 1)— th sheet through the point (A2fc+i, 0) pass anti-clockwise 
to the iV-th sheet through the point (A2fc, 0) and return to the initial point through the iV-th sheet. 
We remark that on Figure 1, when N > 3, the /3-cycles placed from the second to the (N — 2)th 
sheet should intersect the cuts only in the branch points. It is possible to drop this requirement 
and plot appropriate number of loops around each branch point. 

The action of the automorphism J on the basis of cycles is given by 

J(oii+ sm ) = ai + ( s+ i) m , i = l,...,m, s = 0, . . . ,N - 3, 

JV-2 

J(ui+(N-2) m ) = ~ / ] OL i+sm , i = l,...,m, 

s=0 

J(Pi+sm) — A+(s+l)m — A) S = 0, . . . ,N — 3, 

J(fli+(N-2)m) = -A) i = l,...,m. 
The basis of canonical holomorphic differentials reads 

(2.4) du j+sm (P) = - j^ dA, j = l,...,m, fl = 0,..., AT -2. 
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The induced action of J on the holomorphic differential is given by 

J(du j+sm (P)) := du j+sm (J(P)) = AJ . +1 g ff' dA, i = l,...,m, « = 0, . . . , N - 2, 

where p is the N-th root of unity. The (N — \)m x (N — l)m matrices ^4 of a-periods and B of 
/3-periods are expressible in terms of m x m-matrices 



(2.5) A s+ i = i> duk+ms , B s+ i = i> du k + ms , s = 0,...,N-2, 

V a * J k,j = l,...,m ) fe,i=l,...,m 

in the following way. Let us introduce the (TV — l)m x (N — l)m dimensional matrices 

(2.6) n A 



;,fc=i,...,jv-i 



„-i(fc-l) _ -i(jV-X)' 



,fc=l,...,iV-l 



Then 



(2.8) 4 = \<bdu k \ =T>i&g(A l ,A 2 ,...,A N - l )1l A , 




(2.9) B = \<b du k \ = Diag(Bi, B 2 , . . . , B S -i)n Bt 



where Diag(_4 1; _4 2 , ■ ■ ■ ,An-i), Diag(Si,S 2 , . . . ,Bjv_i) are the block diagonal matrices having 
as entries the matrices A s and B s , s = 1, . . . , N — 1, respectively. 

The normalized holomorphic differentials di> = (dwi, . . . , duny— i)m) are given as di> = duA -1 
and the Riemann period matrix II, 



n = I d> dvj 



j,k=l,...,(N-l)m 



is given by II = lZ~ A 1 D'mg(Ai 1 Bi 1 A^ 2 1 B2, ■ ■ ■ jA^^Bn-^TZb with 1Z A and TZb defined in (|2.6|l 
and H2.7|l . The matrix II is necessarily symmetric and has positively defined imaginary part. The 
space of such matrices is called Siegel half-space H. g . Denote also Jacobi variety of the curve 

Jac(C) = c 9 /i 9 ©n. 

To complete this Section we recall the important Ranch variation formulae |33| for the period 
matrix LT 

(2.10) An, = 2^ Res {^E^^)^(^)} - 

where i, j = 1, . . . , 2m, k = 1, . . . , 2m + 1. We remark that in the case of hyperelliptic curves the 
formula l|2.10|l appeared already in the Thomae article For general surfaces the infinitesimal 
variation of Abelian differentials and their periods with respect to Beltrami differentials is due 
to Fay 53 ■ Korotkin |2Zj reduced the variation formula to useful form p.lOjl . The proof of the 
variational formula for the holomorphic differentials can be found in |22j . 
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2.1. Theta functions. Any point e e Jac(C) can be written uniquely as e = (e, 8) (n)> where 

e, 6 e 1R 9 are the characteristics of e. We shall use the notation [e] = [f]. If e and <5 arc half 
integer, then we say that the corresponding characteristics [e] are half-integer. The half-integer 
characteristics are odd or even, whenever A(S, e) is equal to 1 or modulo 2. Here and below 
the brackets ( , ) denote the standard Euclidean scalar product. The Riemann 9- function with 
characteristics [e] is defined in H g x Jac(C) as the Fourier series 

(2.11) 9 [e] (z; II) = exp(m (nil + 511, n + 6) + 2m (z + e,n + 5)) 

= 6{z + e + SH; IT)cxp{2i7r(<5, (z + ^511)) + 2nr(e, 6)}. 

For [e] = [0] we write 6>[0](z; II) = 9(z; II). The 0-function is an entire in z e C 9 and has periodicity 
properties: 

(2.12) 9 [e] (z + m! + mil; II) = cxp{-2z7r(m', (z + ^m'U)) - 2nr(e, m')}0 [e] (z; II), 

where m and m' are integer vectors. 

The ^-function with arbitrary characteristics satisfies the heat equation 

(2.13) ^-^-e[e](z;TL) = 2«r(l + 5 ktl )-=-0[e]{z;n), k,l = l,...,g. 

The zeros of the ^-function are described by the fundamental Riemann singularity theorem. 

Theorem 2.1. Let e G Jac(C) be an arbitrary vector and Q € C- arbitrary point. Then the 
multi-valued function 

p^e^J dv - e; r/j 

has on C exactly g zeros Qi,Q2, ■ ■ ■ ,Qg provided it does not vanish identically and 

(2.14) e = E/ dv - K Qo> 

i=i J Qo 

where Kq q is the vector of Riemann constants 

(2.15) {KQoh = i±Ikl- J2 i Av t (P) f Av . 

Furthermore, the divisor Yli=i Qi * s non special. 

For a point P € C, we define the Abel map 21 : C — ► Jac(C) by setting 

(2.16) 21(f) = / dv, 



JQo 

for some base point Qo € C. For a positive divisor V of degree n the Abel map reads 

r v 

dv. 



21(2?) = f , 



There exists a non-positive divisor A of degree g — 1 such that 
(2.17) fZ(A-(g-l)Q ) = K Qo , 
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where Kq has been defined in Ij2.15|l . The divisor A is called the Riemann divisor and satisfies 
the condition 2 A = ICc, where JCq is the canonical class (that is the class of divisors of Abelian 
differentials). The vector e defined in i|2.14[l can be written in the form 



(2.18) 



e = 2l(^Qi-Q -A). 



i=l 



3. Characteristics supported on branch points 

In this section we are going to compute the characteristics [14k] of the Abelian images of branch 
points P k = (A fe , 0) 



u k = J dv, k = l,. 

in terms of the period matrix II. 

Lemma 3.1. The relations 

rPzk-i N — I - s 

dVk+sm — 



, 2m + 1, Poo = (oo,oo), 



P2k + 1 

Plk + 2 
P2j 

P2 3 +l 



N 



P2 k+ i N-l-s 
dwfc+sm — - 

-P2fc + 2 



N 



dv 



dvk+ 



j+sm =0, j ^ k,k+ 1, j = 1, . 

N-2 

N 



, m, 



N-l. 



N 



-n 



k-\-sm,j 



N-2 



k+smj+rm 



r=l 



are valid for k. j — 1, . . . , m, s = 0, . . . , N — 2. 



The proof represents a generalization to the case N > 2, of the derivation of half- integer char- 
acteristics given e.g. in |16|: for detail see the proof of this lemma in jllj . 

From the relations given in the Lemma |3. II we are able to write the characteristics [14 k] in the 
form 



2m+l\ — 



ml 



AT 




— 



(N-l)ml 



N-l 
N 







[U 2 



m I 
1 



N 



sml 

_J_ 

N 
s 

N 



(JV-l)mJ, 
"AT 



u N 



k I 


A' 



J_ 

N 



1_ 

N 



fc+(s-l)m| 



JV 



N 



J_ 

N 
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k+(N-2)ml 



JV-1 
N 



J_ 

N 



J_ 

N 



[U 



2k 



fej 

_J_ 

N 

_L 

N 



N 



o 



J_ 

N 



fc+(a-l)mj. 

_ J_ _J_ 
N N 





J_ 

N 



k+{N-2)ml 
1 



N 



J_ 

N 



J_ 

N 



[U 2 



f 

N 



N 



JV 



l+(s-l)ml 

_ J_ _J_ 

N N 





N 



_ J_ 

iV-1 
N 



J_ 

N 



J_ 

N 



[Ui 



N N 





N 



N N 





N 



N N 





JV 



In the following wc determine the vector of Ricmann constants and Riemann divisor of the curve 

C. 

Lemma 3.2. The vector of Riemann constants computed in the homology basis described in Fig- 
ure 1 and with base point P^ eguals 

(3.1) Koo = (JV - 1) J2 / dv > P 2fe = (A 2fc ,0). 

fc=ip 

■IOC 

The Riemann divisor A of the curve C in the homology basis described in Figure^is equivalent to 

m 

(3.2) A = (7V-l)^P 2fe -P oc . 

fe=i 

Proof. The proof of (|3.1fl is obtained by direct calculations from the definition (|2.15|) and Lemma l3~Tl 
The relation l|3.2(l follows immediately from l|3.1|l . □ 

Following Diez [7|, we describe a family of non-special divisors of degree g on C, supported on 
the branch points. Denote by s = (si, . . . , S2 m +i) a 2m + 1-vector with non-negative entries s&, 
satisfying the condition 

2m+l 

(3.3) ^ Sj = (JV — l)m, < < JV - 1. 

i=l 

To each vector s we put into correspondence the divisors 

2m+l 

(3.4) V s = s kPk 

fc=i 



s 
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where Pk = (Xk, 0), k = 1 . . . , 2m + 1, are branch points. In particular we shall consider the divisor 

— ,„ t (^)^ 

(3.5) V m = (N - l)P n +... + {N- l)P im _ 1 + (N - l)P im , 

where the indices «2, ■ ■ ■ , 4} C {1, . . . , 2m + 2}. Among the divisors with m+ 1 branch points 
we consider the divisor class D m |i which contains |(m + 2)(m + 1) I ^ ^ J divisors 

(3.6) = (N — l)P h +... + (N- l)P im _ 1 + (JV - 2)P im + P Wl , 

where the indices «2, ■ ■ • , *m+i} C {1, . . . , 2m + 2}. It is out of the scope of the present 
manuscript to classify all the non-special divisors of the form (|3.4|) . However we can single out two 
families of non-special divisors. 

Lemma 3.3. The divisors T> m defined in |5*. 5)) are non-special and the divisors P m ^i defined in 
\3.b\l are non-special for N > 3. At N — 3 the divisors 

(3-7) V m+1 =2P tl +... + 2P lm _ x + P im + P im+1 , 

are non-special if one of the following conditions are satisfied 

• i m and i m +i have different parity, 

• i m and i k have different parity for k < m, m > 1, 

• ik and ij have different parity for k,j<m and m > 1. 

Proof. Assume the opposite: suppose that the divisor T> m or T> m _^_i is special, this means that 
there exists a non-constant meromorphic function f(X,n) whose divisor of poles is T> m or X , rrl _|_ 1 . 
For simplify the proof, we assume that Pi k ^ P^, k = l,...,m + l. Let I; be a subset of 
{1, 3,5,..., 2m + 1} with I distinct elements. The function 

<P(A,y) = f(X,y) n( A - A 0> l = m,m + l, 

ij Gli 

has poles only at infinity. It follows from the Weierstrass gap theorem, that the ring of meromorphic 
functions with poles at infinity is generated in the case of the curve y N — p(X)q N ^ 1 (X) by powers 
of A and functions y k / 'q(X) k ^ 1 , k = 1, . . . ,N— 1 with q(X) = rij=i( A ~ A 2j)- Therefore the function 
4>{X, y) can be written in the form 

N-i k 

(3-8) 0(A,y)-P o (A)+^P fe (A) n ^ 7 -, 

k=l Q ^ ' 



where Rk(X) are polynomials in A. 1 



We remark that ordoo ^Pfc(A) k _ 1 J ^= ordoo lPj(A) ^ or ^ ^ because otherwise 

(3.9) A^ord A P fc (A) + k = Nord\Rj(X) + j, 

which implies k — j. From this observation it follows that 

(3-10) ordoo (/(A, y) 1] ( A ~ A *J) = ord - MU A )^T7 



gi-i(A) 



4n this point our proof differs from that given in |7| which is working for Galois covers of the form y 1 * 
rifcLi(^ ~ ^k) where the ansatz for the function 13.81 can be written as RiV 1 ■ 
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for some < j < N — 1 . Moreover 

ord^ifi^y) JJ(A-A <B )) = -JV|Ii| + fcj J l = m,m+l, 

in&l 

where ki, I = m^m + 1, is the order at infinity of /(A, y) and = I. From the equation of the 
curve we get deg y = mN + 1 . Therefore the equality 1)3.10(1 can be written as 

JV|Iz| - ki = N(r 3 +m)+j, I = to, m + 1, 

where rj is the degree of Rj(X), so that 

j = N(\Ii\-r j -m)-k l > 0. 

When I = m that is = to, it follow that rj = 0, j = 0, k m = and 

/(A ' v) = n<„ e jA-Ao 

and contradicts the assumption that /(A, y) has divisor T> m . 

When I = to + 1, that is |Ij| = m + 1, two possibilities occurs: (i) rj =0, j = N — k m+ i, < 
fcm+i < N and (ii) rj = 1, ki = 0, j = 0. This latter case can be easily excluded while for the 
former one we have 

I y N-k m+1 

/(A,yH n, j£ i m+1 (A-A. tj ) g ^ +1 -HA) ' ord °°(/(A,2/)) = W, 
which has divisor 

ro+l m 

Div/(A, y ) = -AT ^ P in +{N- k m+1 ) J2 p 2j+i + fcm+i I] 

Namely the divisors of poles of /(A, y) is 

Div po i os /(A, y) = (N — k m+ i) ^ Pi n +k m +i ^2 Pi *> 

in&m+i^n even i n ,2„, odd 

and for N > 3, differs from P m +i. This contradicts the assumption unless / is constant. For 
N = 3 the divisor of poles of /(A, y) coincides with D m |i in the following two cases: 

m— 1 

X> m+ i = 2 P lfc + P im + P im+1 , 

k=l 

with 

i m , im+i € {2, 4, 6 ... , 2ra}, i fc € {1, 3, 5, . . . , 2to + 1}, k = 1, . . . , to - 1 

or 

4, Wi G {1, 3, 5 . . . , 2to + 1}, i fc e {2, 4, 6, . . . , 2m}, k = 1, . . . , to - 1. 
We conclude that the divisors (|3.7I) are non special when one of the following conditions are 
satisfied: 1) i m and i m +i have different parity; 2) i m and ik, k < m, m > 1 have different parity; 
3) and have different parity for j, k < to, to > 1. □ 

By the above lemma and by the Riemann singularity theorem l2.ll the vector 

e m = m(V m - Q - A), Q (£V m . 

is a non singular 1/N period, namely 9(e m ]lV) ^ 0. In the above expression we can get rid of the 
base point Qq. For the purpose, we introduce the divisor 

V = Q + J(Q) + --- + J N ~ 1 Q, 
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which is independent from the point Q and satisfies the relation 

V = NPj, j = 1,..., 2m + 2, P 2m+2 = -Poo- 
Let Jo = {2, 4, . . . , 2m + 2} and Iq = {1,3,..., 2m + 1} be a partition of the branch points in odd 
and even and P^ = Pi m +i- Let Ji C Jo and 1\ C lo be a partition such that 

|Ji| + |/i|=m + l. 

Consider the vector 

(3.11) e m = W((N - 1) Pi + (N - 1) Yl P J ~ V ~ A )- 

Then from Lemma 13.31 the vector e m is non singular, namely 9(e m ;H) ^ 0. In the same way we 
define the vector obtained from the divisors T> m +i in the following way 

(3.12) e m+1 = a((JV - 1) J2 P i + ( N ~ !) P i + ( N ~ 2 ) P< ™ + P W» _ V ~ A )' 

ie/i jGJi 

where now 

|J X | + |Ji| =m- 1, JiC/o, Ji C J , 

and 

j m e (Jo - h), 4+i e (Jo - Ji)- 

From Lemma 13.31 the vector e m ^i is non singular. 

In the next sections, we are going to study the Szego kernel associated to the characteristics 
jSZEHl and (332)1 . 

4. Kernel-forms 

The Schottky-Klein prime form E(P, Q), P, Q € C is a skew-symmetric (— 4, — ij-form on C x C 

0[ 7 ] (fdv-M 

(4.1) £(P,Q) = 



E?=i 4^[7](o ; n)d^(P)^Ef=i £-ehmu)dv 3 (Q) 

where [7] are non-singular odd half-integer characteristics. 

The prime form does not depend on the characteristics [7]. The automorphic factors of the 
prime form along all the cycles ak are trivial; the automorphic factor along each (3k cycle in the Q 
variable equals exp{— irillkk — 2ixi J p dt^}. If the points P and Q are placed in the vicinity of the 
point R with local coordinate x, x(R) = 0, then the prime form has the following local behavior 
as Q -> P 

(4.2) E{PQ) = ^h X } P) (1 + 0(1))). 

The canonical bimeromorphic differential uj(P, Q) is defined as a symmetric 2-differential, 

(4.3) u(P,Q) =d x{P) d xm log E(P,Q). 

All a-periods of uj(P,Q) with respect to any of its two variables vanish. The period of the 2- 
differential oj(P, Q) with respect to the variable P or Q, along the f3k cycle, is equal to 2mAvk{Q) 
or 2indvk{P) respectively. 
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The 2-differential cj(P, Q) has a double pole along the diagonal with the following local behavior 

[HI E] 

(4.4) lo{P, Q) = ( -. 1 + H(x(P), x(Q)) + higher order terms ] dx(P)dx(Q), 

where H(x(P), x(Q))dx(P)dx(Q) is the non-singular part of u>(P, Q) in each coordinate chart. The 
restriction of H on the diagonal is the Bergman projective connection (see for example |37j ) 

(4.5) R{x(P)) = 6H(x(P),x(P)), 

which depends non-trivially on the chosen system of local coordinate x(P). Namely the projective 
connection transforms as follows with respect to a change of local coordinates x — > f(x) 

R{x) -> R(f(x))[f(x)] 2 +{f(x),x}, 

where { , } is the Schwarzian derivative. 

The Szego kernel S [e] (P, Q) is defined for all non-singular characteristics [e] as the (^, i)-form 
on C x C Q3] 

0[e] ljdv,u\ 

(4 - 6) s ^ p ^ = mmmQy 

The Szego kernel transforms when the variable Q goes around a k and /3fc-cycles as follows 

S[e](P, Q + a k ) = e 2 ^S[e](P,Q), 
(4.7) „ k = l,...,g 

S [e] (P,Q + f3 k ) = e~ 2m,ik S[e} {P, Q), 

The local behaviour of the Szego kernel when Q — > P is 



(4.8) 5 [e] (P, Q) = v|^g^ [1 + T(o ; (Q))(x(Q) - x(P)) + 0((x(Q) - x(P) 2 )] , 
where 

(4.9) T(x(Q))dx(Q) =J2q^ log0[e](O;n)d Ufc (a;(Q)). 

k=l Zk 

The important relation |14| . Cor. 2.12, connects the Szego kernel with canonical bimeromorphic 
differential 

(4.10) S[e](P,Q)S[-e](P,Q)=cj(P,Q)+ £ logtf [e] (0; n)dv k (P)d Vl (Q). 

/c,/— 1 

In the following we are going to give an algebraic expression for the Szego kernel associated with 
the characteristics 

e m = W((N - 1) Pi + (N - 1) Yl P o ~ V ~ A )' 

ieh jeJi 

defined in (|3.11|) . For simplicity, we assume 2m + 2 ^ J\ . We define the functions 

(4.11) MP,Q) = X f2~ X * , k = l,...,2m+l. 

x(P) - Afe 
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Theorem 4.1. The Szego kernel with characteristics [e m ] defined in 13.11$ is given by the formula 



S[e m }(P,Q) 



(4.12) 



y/dx(P)dx(Q) 
N(x(Q) - x(P)) 

N-1 i o 



E 



N-1 _ 



rw 2 $i 



y-i r £+1 1 ■ 

"373 1 ~TT J 



where { } is the fractional part and 

J 2 = Jq — J\ 



2m + 2, I 2 =I -h. 



In particular, the Szego kernel with zero characteristics is obtained by fixing I\ — Iq and J\ 
and takes the form 



(4.13) 



S[0](P,Q) 



1 v /dx(P)dx(Q) 



N-1 



N x(Q)-x(P) ^ 



^ ( q{x{P))p{x{Q)) \ 



,p(x(P)) q(x(Q))J 
where the polynomials p(X) and q(X) have been defined in 



The proof is based on the uniqueness of the Szego kernel (see e.g. Narasimhan [32]) and the 
results of Nakayashiki 31 . Indeed it is sufficient to check that the expression H4.12[l is regular 
everywhere but on the diagonal where P — Q and that its divisor in the variables P and Q 
coincides with the divisor of the Szego kernel given by the formula l|4. fijl . Using the coordinate 
chart given in (|2.3|) , the regularity of the expression l|4.12|l can be checked in a straighforward way. 
In the same way, by first fixing P — P^ and then Q = P^ with ii 6 Jo U Jo it is straighforward to 
obtain the divisor class of ^4.12|l . 

Corollary 4.2. The expansion of the Szego kernel \4-l%fy along the diagonal takes the form 
S[e m ](P,Q) 



y/dx(P)dx(Q) 

x{Q) - x{P) 



1 



I.+ I ^MQ),Q} + ^e m ](Q) ) (x(P)-x(Q)) 2 



where the function <p[e 
<f>[e m ](X) = 

(4.14) 



m ] (A) takes the form 
N 2 — 1 




Xi) 



d_ rwA 

^ og rw 2 (A-A,) 



rfA ° g IWA-A,) 



5) 



247V 2 



rfA ^n^CA-A,) 



^A ^rWA-A,) 



The proof of the above corollary is obtained by direct calculation. We remark that for Ii 
and J\ = the second term in l|4.14|) disappears and we obtain 

2 



flop) = 



N 



1 



247V 2 

where p(X) and q(X) have been defined in l|2.2|l . 
Since 

-e m = 2t((iV - 1) J2 P >- 

ieio-ii 



dX ° g q{X) 



(N-1) J2 Pj-V-k) 
jeJo-Ji 
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because 

(N - 1) p k - 2X> - 2A = 

we conclude that 

4>[e m ]{Q) = 4>[-e m ](Q). 
Therefore, the following identity holds 

g [eJ(P,Q)5[-e m ](P,Q) = ( ^ ) P )^ )2 



(x(P) - x(Q)) 2 + 



x |l+ ~{x(P),P} + 2<p[e m ](P) 
The relations (|4.8(l . Ij4.9jl and the corollary PT2l imply 

(4.15) ^-e[e m ](^n)U =0 = 0, fc = l,..., 5 , 

for the vector e m defined in 

In the following we are going to give an algebraic expression for the Szego kernel associated to 
the characteristics 

(4.16) e m+1 = 2l((iV - 1) Y, p i + ( N ~ 1) £ P i + ( N ~ 2 ^ + P <« ~ V ~ A )' 

ie/i jeJi 

where 

| Ji| + | Ji| = m - 1, i m G J - Ji, jm e Jo - Ji 

The construction is very similar to the previous case. Also in this case, for simplifying the notation, 
we assume that 2m + 2 ^ J\ and j m ^ 2m + 2. We define the following sets 

h = Jo - Ji - i m , <h = Jo - J\ - jm -2m + 2. 

Theorem 4.3. The Szego kernel with characteristics [e m _|_i] defined in 14-16\ ) is given by the 
formula 

/W](P,Q) = v^)^) . 



N(x(Q)-x(P)) 



(4-17) 

E 



s=0 



n,\2N iv/TT /\2JV { it } / , \ — J 



where { } is the fractional part and the function ipi has been defined in \4- 1 

The proof is based on the uniqueness of the Szego kernel (see e.g. Narasimhan [221 ) an d the 
results of Nakayashiki [2] . 

Corollary 4.4. The expansion of the Szego kernel along the diagonal takes the form 



S[e m+1 ](P,Q) - x{Q) _ x{p) 

x [l + (^{x(Q),Q} + <P[e m+1 ](Q) ) (x(P) - x(Q)) 2 
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where the function <p[e m +i]{X) is 

N 2 - 1 

0[e™+i](A) 



(4.18) 



247V 2 

d_, . (A-A im ) 
dX ° g (A - Xj m ) 



d_, rWA-A^) 

rfA ^iwa-a,) 



2" 



IWA-AQ 

rfA ° s n ieJl (A-A,) 



2(AT- l)(iV-5) 
24iV 2 

2(iV 2 - 12iV + 23) 
24iV 2 



rwA-^) 
rfA og n JGJ2 (A-A 3 ) 



rfA og n i6Jl (A-A 3 o 



d, (A — Ai ) 
log 



dA 



(A - Aj m ) 



where 



I{ = h + 4, J' 2 = -h + j r , 



The proof of the above corollary is obtained by direct calculation. 
Since 

-e m+1 = 2l(iV - 1) J] P + (JV - 1) ^ Pj + (N- 2)P lr 

i€.h j£Jo—Ji-jm 

because 

(N - 1) ^ P fc - 2X> - 2A = 

fceiouJo 

we conclude that 

4>[e m+x ]{Q) = (j>[-e m+1 ](Q). 
Therefore, the following identity holds 

dx(P)dx(Q) 



P jm - V - A) 



S[e m+1 ](P,Q)S[-e m+1 ]{P,Q) 



x 1 



-{x(P),P} + 2#^ +1 ](P) 
o 

We remark that from 14. HI . (|4.9|l and the corollary l4"2l 



(x(P) - x(0)) 2 
(z(P)-z(Q)) 2 



(4.19) 



^-0[e m+ i](z;II)| z= o = O, k = 1, 



for the characteristics e m +\ defined in (|3.12l) . 

5. Algebraic realization of the canoncal bimeromorphic differential 

The canoncal bimeromorphic differential can be given in an algebraic form due to Klein |23 | .|24 | 
also an d H31- To develop this approach we first write the third kind differential with poles 

in two arbitrary points. For the purpose we consider an arbitrary curve C given by the polynomial 
equation /(A, fi) = of degree N in the variable /i. We suppose that the curve C has a branch 
point at infinity. Let 

(5.1) *(A,/i) = (l,Vi(A,/i),...,Viv-i(A, ) t 1 )) 

be the basis in the ring 0(C) of meromorphic functions on C with the only pole at infinity. There 
exists a vector function 



(5.2) 



#(A,/i) = (l,$i(A,/i), 



(A,M)) 
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for which 

(5.3) (*(*,„), *(A,„')) = fiX ^']-_ fiX ^\ 

Let Q — (A',//) and R = (A",//') be two arbitrary points of the curve C. The third kind 
differential Q,qji(P) with simple poles with residues ±1 in the points P = Q and P — R is given 
by the formula 

(5 ' 4) " Q ^ (P) - {—X^\i A^A^J j^)' 

with / M (A,/x) = df(X,fi)/dfi. In the case of the curve (12.1(1 the vectors \1>(A, /x) (see the proof of 
Lemma \'S. 31 ) and 3?(A, /i) take the form 

(5.5) *:A./n (l,fi 



g(A)'"-'g(A)^ 

(5.6) #(A, ^) = (n N -\ M^" 2 , aCA)^- 3 , • ■ • , g(A)"- 2 ) . 

The canoncal bimeromorphic differential can be obtained by differentiation of l|5.4|l as follows 

„ (P , c) . d ^fflM_^_ +dx(P , c) , 

where d%(P, Q) is the 2-form uniquely defined by the requirement that 1)5.7(1 is a symmetric bi- 
diffcrential normalized with respect to the a cycles and with the only pole of second order along 
the diagonal. The above arguments lead to the following algebraic expression for the canoncal 
bimeromorphic differential: 



(5.8) W (P,Q) ° ' 



dX' A - A' 



y Vg(A)- 1 
^ /x s g(A') s_1 



dAdA' , /n „. 
+d x (P,Q), 



(5.9) d X (P,Q) = ~EE^^t^(W 

S = l j = l ^ ^ 

with certain polynomials lZ s _j(X, /i). Here we do not need the exact form of these polynomials. 

Expanding the above expression along the diagonal, we obtain the Bergman projective connec- 
tion 

«(P)=i W P),p }+ IS$ 



4 JV ^ p(A) g(A) y 127V 2 \dX b g(A) 

Combining the above relation with the expansion 14.1511 of the Szego kernel along the diagonal and 
the Fay relation ((4.1U|) which connects the canoncal bimeromorphic differential and Szego kernel, 
we obtain an algebraic expression for the second derivatives of the theta function, namely 

1 9 d 2 

7r7r log^[e m ](0;n)d Ufe (P)d« / (P) 



(dx(P)) 2 ^ dZkdZi 

(5.11) 



IN-lf^p(X) feg(A)\ N 2 -l fd ,„p(A)^ 



A N \ p(X) q(X) J 127V 2 \dX b q(X) 

where <f)[e m ](P) has been defined in 1(4. 14|) . 



iog^ +20[e m ](P), 
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The following relation will be useful for the proof of the Thomae type formula. It connects 
the function dx(P, Q) defined in l|5.9|) with the derivative with respect to the branch points of the 
determinant of the matrix A of a-periods. 

Lemma 5.1. For s = 1, . . . , N — 1 the following identities are valid 

d 1 m 
(5.12) logdet^^^j »' i = l,...,2m+l. 

II (A, - Xi) J' =1 

Proof. We integrate w(P, Q) in the variable P along the a/c-cycle and expand § w(P, Q) = in 
the variable Q in the vicinity of the branch point (A^ , 0) where the local coordinate is introduced 
as x(Q) = Xi + £ , i € {1, . . . , 2m +1}. In this way we obtain for every fixed s the rule of variation 
with respect to the branch points of the first line of .A s -block of the period matrix A: 

d 1 ™ 
Tf-(A s )i,k = -j^i y2 R -.- :A - k=l,...,m. 

The derivative of all other lines can be obtained from the equivalence 

n— l 



^(A s ) n , k = Xr 1 ^M s ) hk + ±Y,( A °hkXr l -\ n = 2 ) ...,m, fc=l,...,i 
Therefore one can write 



o m 

_(A) = V s (Xi)A s , TrP = £ A r 1 ^ J (A l ). 
' j'=i 

The equality (|5T2l follows. □ 
Proposition 5.2. TTie following relation is valid 

(5.13) A logdet.4 = - Res ?* ( f' ^ ■ 

dXi P=Pi (da;(P)) 2 

Proof. The residue of l|5.9ll is 

1 N-l m 

to EEr^iiA,). 

n (Ai - ao -=i j=i 

Then use the decomposition (|2.8|) and the Lemma IS~Tl □ 
Example 5.3. In the case of elliptic functions the formula (|5.13|l represents the known relation 

du_ _ 1 1] + XjLO _ 

dx t - 2{x i -x j ){x i -x k y '^^-m,* 

where the Weierstrass notations are used. 



THOMAE FORMULAE FOR SINGULAR Z N CURVES 



17 



6. Derivation of the Thomae type formula 
Now we are in a position to derive Thomae type formula 
Theorem 6.1. Let [e m ] be nonsingular 1/N characteristics, given by the formula 

e m = St((JV - 1) p i + ( N ~ !) P o ~ V ~ A )< 

ieli jeJi 

where J\ C Jo = {2, 4, . . . , 2m + 2} and I\ C Iq = {1, 3, . . . , 2m + 1}, with 

|Ji| + |/i|=m + l. 
The Thomae type formula takes the form 
(6.1) 

0[e m ](O;II) 



AN lli=l at - z -^i TT f\ \ \N(N-1) TT (\ a \N(N-1) 

~ ( 2m )2mN(N-l) 11 (M-Mk) [I IAm+1 - A 2fe+ i) 

V 7 Kk 



i<k<m 0<i<fc<m 

2(Af-l) 



X 



ii- / ,, ./ ' :A - - a j) n ie / 2je j 2 (^ - a j 



,Ilie/i,ifcei2^ _ A fe) n je j 1: fcej 2 ( A j _ A fe) y 
where 

h = Io-h, -h = Jo - J\ - 2m + 2. 
When I\ — Iq and J\ = 0, i/ie above formula reduces to 

(6.2) 0[o](O;nr = |g^g n (^"W^" 1 ' n (W-W)^. 

^ 7 l<j<fc<m 0<i<fe<m 

Proof. Using the heat equation (|2.13[) . the Rauch variational formula (|2.1(J|) and (|4.15|l we have 



(AT-I)jti 

(6.3) ^\ og e[e m ]{Q-Ii) = Y, ^log^[e m ](0;n) 



(JV-l)m 2 AT 



o\i oIL kr dX 

K,r— 1 



= 9 D ^f ^TTTmW E ^-log^[e m ](0;n)^d^(pW)d« r (pW) 
2 p=(A 4 ,o) (dx(P)y ' dz k dz r ^— ' 

To proceed we shall represent the residue as a logarithmic derivative. The equality l|5.11|l and 
the proposition 15 . 21 enable us to compute the residue in the r.h.s of l|6.3|l obtaining 

^-log0[e m ](O;II) 

oXi 

8 a=a, \ p(X) q(X) I 2 aXi 

+ p Res o) ^e m ] (P) - ^ log ^|) 2 , 

where <p[e m ](P) is defined in (|4. 14|l . 

Let us compute each residue in the above formula: 

Re S ^\ = 2— log TT (X 2l - X 2k ) and Res ^ = 0, 

A=A 2l p(X) 3X 2i 11 A=A 2l + 1 p(X) 

1 < i < k < m 
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Res / \ \ — = 2 — log 



A=A 2 i+i <?(A) 



<9A 



2i+l 



J| (A 2i +i - A 2 fe+i) and Res aA J^ 

1 < i < k < m 



Also we have for any partition I U J and any A 

2 



d n(A-A,) 



a n (Ai-Afe) n (Ai-Ak) 

° S II (Ai-A,-) ' 



Furthermore 



Res 

A=A„ 



d_ ] n.e/^A-A,) 
dX ° S n ie , 2 (A-A,) / 



A log n^(A-A.) 

rfA g IWA-A,) 



___ 

<9A„ 



log 



n te / 1: fcgj 2 (^ ~ Afc)IIjgj 1 fcgj 2 (Aj - Afc) 
riie/i je./i(Ai ^ a j) n, e i 2 ,jej 2 (Ai - Aj) 



N(N-l) 



The integration in A„ then gives 
(6.4) 

% m ](o ; n) 4Ar =c 2Ar n detA 2jv n (** - A 2fc ) Ar(Ar - i) n (w-a 2 , + i) 

2=1 l<z'<fc<m 0<?<fc<m 

/ Ilig/u-gJxCAi - -V. ! 11,. ;, ., ,;/A, - Aj) \ 
\n i g/ 1 ,fcg/ 2 (A i - Afe)IIiGJi,feej 2 (Aj - Afc) y 

where C is a constant independent of II and the branch points. To derive the above relation we 
also used the decomposition of the period matrix A into blocks A% ■ 

To compute C we shall use the original Thomae arguments [HEj • We consider the Thomae type 
formula (|6.4(1 at zero characteristics, that is, the partition I\ = Iq and J\ = and pinch the branch 
points in the following way 



A 2 fc = Vk 



A2fc-i=«fc + e k=l,...,m, < e < 1. 



In this case the l.h.s of (|6.4|l becomes 0[O](O;II) = 1 + 0(e). Regarding the r.h.s the following 
relations are needed: 

j-i 

lim(A)i,j = 2m— -i 



II ( v i ~Vk)(Vi - A 2m+ i)« 

k=ii 
fe=l 



so that 
(6.5) 



limfdet A s ) = (2mY 



1 



1 



f] (Vk ~Vj) J] (v k - V2m+l)~ 
k<j k=l 
k,j=l 



Combining the above relations the expression for the constant C = (2m) 



— m(JV-l) 



□ 



It is straightforward to check that for N = 2 the formula (|6.1(l coincides with the original 
Thomae formula 1)1.1(1 . 

To complete the section we give the Thomae formula for the characteristics [e TO +i] defined in 

irm . 
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Theorem 6.2. Let [e m ^i] are nonsingular 1/N characteristics, given by the formula 

e m+ i = 2l((iV + P i +( N ~ 2 ) P ^ + p i m — D — A), 

ie/i jeJi 

where Ji c Jo = {2, 4, . . . , 2m + 2} - {2m + 2} and h C I = {1, 3, . . . , 2m + 1}, wii/i 

|Ji| + |/i|=m-l 

and j' m G Jo — Ji — 2m + 2, i m 6 Jo — -^l • J^-e Thomae type formula takes the form 
(6.6) 



^ 27r )2mJV(iV-l) 

V ' l<i<fc<7Tl 

Vr[i e j(,fe e / 2 (Ai - Afc)rijeJi,*e^(Aj ~ Afc) y 

/ n je .7 2 (A, m ~A J )n. tg/1 (A t -A J j \ 4(JV ~ 2) 



0<i<k<m 
\ 2(JV-1) 



where 



I[=h+hn, h=h-h-imi J 2 = Jo ~ Ji - j m - 2m + 2, Jj =J 2 +.?',„• 
The proof of the above formula follows the steps of the proof of 1)6. 1|) . 

Example 6.3. We consider as an example the trigonal curve of genus two (studied by Hutchinson 
|19| . see also \TI\ for more details) 

(6.7) A1 3 = (A-A 1 )(A-A 3 )(A-A 2 ) 2 . 

In the homology basis given in Figure ^ the Ricmann period matrix is of the form 



(6.8) 



n = 



2T T 
T 2T 



T = 



n/3F(± f;l;l-f) 
3 F(\,\;l\t) 



where F(a, b, c; t) is standard hypergeometric function and t = (X% — Ai)/(A3 — Ai). In the homology 
basis given in figure Q the characteristics of the branch points are 



[2ti] = 



2 2 

3 3 




ra 2 



2 2 

I I 

3 3 



[a 3 



o 

1 2 

3 3 J 



[2U 



[0], 



whilst characteristics of the vector of Riemann constants are [K p^] = 2[2l2]- Given 

7 = {1,3}, Jo = {2,4}, P 4 = Poo, 
we introduce the following partitions of the branch points. 

(1) The partition 

h ={1,3}, Jx = or h = 0, J a = {2,4}, 

which corresponds to the characteristics [e m ] = [0], 

(2) The partition 

7i={l}, Ji = {2}, 

or 

Ji =0, Ji = 0, i m = 3, j m = 2, 



which corresponds to the characteristics [e m ] = 



l l 

3 3 
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(3) The partition 



Ji = {3}, Ji = 2, 



or 



Ji = 0, Ji = 0, i m = l, j 
which corresponds to the characteristics [e m ] = 
In the first case, the Thomac formula takes the form 





2 1 
L 3 3 



where 
(6.9) 



3 (As-Ai) 1 ^ ^3'3 



A 2 



where p 3 — 1. The second partition gives 



and the third partition 



l l 

3 3 







2 1 

3 3 



(o ; n) 12 = 



(o ; n) 



12 



AiA-2 
4tt 2 



MM 

At: 2 



(Ai - A 2 ) 4 (Ai - A 3 ) 2 



(A 3 -A 2 ) 4 (A 3 -Ai) 2 



The quotient of the first and third relations leads to 



ft 

V 





2 1 

3 3 



0(0; n) 



(o ; n)^ 



(A 3 - A 2 ) 2 
(Aa-Ai)2 



which we mentioned in 



7. Conclusion 



In this paper we presented the derivation of the Thomae type formula for singular curves with Zpj 
symmetry. We computed explicitly the Riemann period matrix of the curve in the fixed homology 
basis given in figure ^ taking into account the action of automorphism. We also described 1/7V 
periods of the curve in terms of its characteristics given in the form of g x 2-matrices with rational 
entries. 

We considered in the papers a family of non-special divisors supported on the branch points 
and derived the Thomae formula for this family only. But our derivation is general and is working 
for any other family of non-singular 1/A^-periods. The proof given above goes up to the original 
Thomae proof and involved a number of steps such as Rauch variational formulae, calculation of 
the holomorphic projective connection, certain results about the Bergman kernel. But the key- 
point of the proof is the derivation of an algebraic expression for Szego kernel associated to the 
aforementioned 1/N-period. We show that the values of the exponents in the Thomae formula 
follows directly from this expression. 

The results of the paper can be generalized to the family of curves 

H N = (A- Ai) mi ...(A- X k ) m ", mx + ... + m k divisible by N. 

However even for the curve considered a number of open problems remain, among which the 
complete classification of non-special divisors supported on the branch points and the derivation 
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of an algebraic expressions for the associated Szego kernels. We remark that a family which is 
different from the one discussed in the present manuscript has already been studied in jH|. 

Another set of interesting problems consist of the derivation of Thomae formulae for ^-derivatives, 
which generalize the formulae given by Thomae in the case of hyperelliptic curve |36j . Such for- 
mulae are important for obtaining the expressions of the "winding vectors" in terms of ^-constants 
and which are generalization of the Rosenhain formulae known for genus two. 

We can also point out the problem of the derivation for the given curve of the Jacobi-Riemann 
derivative formula, one of the most mysterious ^-formula, which Riemann called as a "pearl of 
mathematics" . 

The answers to these questions are of interest for the applied problems pointed out in the 
Introduction. The investigation of the particular curve undertaken above was stimulated by its 
links with the solvable Riemann-Hilbert problem of rank N considered in 
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